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TARGET TYPE ESTIMATION IN TARGET TRACKING 

Technical field of the invention 

The present invention relates in general to tracking of targets by means of 
measurements from various sensors and in particular to target type estimation by 
using discrete, target type related information. 

Prior art 

In target tracking, data from sensors are used to determine a target track. This 
target track, or entities involved in the process of creating it, is associated with a 
certain target. It is useful if this target could be identified for supporting the 
classification, e.g. friendly aircraft against hostile ones, helping the operator of the 
tracking system to initiate relevant measures. The target type identification may use 
many types of different information, for instance, discrete information associated 
with e.g. ESM (Electronic Support Measures) data, IRST (Infra-Red Search and 
Track) measurements and direct target type observations. There exist many 
different types of targets that may need to be classified, which requires a method 
using efficient algorithms to reduce the computational needs. 

Many multi target tracking systems of known art utilise algorithms which are based 
on a probabilistic approach, why it would be advantageous if extra functionality such 
as target type estimation would retain a probabilistic framework. 

It is, e.g. from S.S. Blackman "Multiple-Target Tracking with Radar Applications", 
Artech House, MA, USA 1986, p. 368-380 known in prior art to use a BayesiaR 
probabilistic framework for estimation of discrete quantities. Target type probabilities 
are suitable estimates, which are easily integrated within the probabilistic framework 
of the total multi target tracking system. Such methods have the disadvantages of 
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requiring more and more data during the calculations when the measurement history 
grows, and this amount of data will eventually prohibit a real time treatment of the 
problem. Another disadvantage with the method is that it exhibits an inherent 
difficulty in treating ignorance or uncertainty. 

It is further known to use e.g. the Dempsters-Shafer method (S.S. Blackman 
"Multiple-Target Tracking with Radar Applications", Artech House, MA, USA 1986, p. 
380-391) in similar applications. One of the disadvantages with such methods is that 
the true probabilistic interpretation is lost. As a consequence, the integration of this 
type of method in a target tracking system of probabilistic type becomes 
troublesome. 

Other approaches, such as neural networks and fuzzy logic are known in the art. 
However, these methods are not easily integrated with other parts of the 
probabilistically based tracking system. 

In the US patent 5,392,225 a method and apparatus for correlating target data is 
disclosed. This method is time static, i.e. uses only information from the last 
measurements. The method computes the likelihood for measurements from 
different sensors to originate from the same target. Disadvantages with this method 
is that neither filtering over time, nor correlation to tracks, is used, which significantly 
limits the degree of accuracy. 

In the US patent 5,392,050 a method of recognising the type of a radar target object 
is disclosed. This method uses a time-frequency analysis of the RCS (Radar Cross 
Section). The method is not applicable to e.g. passive sensors and discrete 
information sources. 

In the US patent 5,282,013 a passive ranging technique for infrared search and 
track systems is disclosed, in which a library of emission vs. contrast is available. By 



WO 98/58274 



PCT/SE98/00956 



3 

measuring the atmospheric background, a type match may be found, which 
eventually gives the target distance. Earlier distance estimates are also used in the 
calculations. However, this* system is integrated in a special type of target tracking 
system and may not be used in a general manner. 

Summary of the invention 

An object of the present invention is to provide a method for target type estimation 
for target tracking purposes using discrete information, which method is possible to 
execute in real time and that is able to handle ambiguities and uncertainties. 

Another object of the present invention is to provide a method for target type 
estimation using discrete information from different types of measurements. 

Yet another object of the present invention is to provide a method for target type 
estimation which easily is integratable in a probabilistic framework for various target 
tracking purposes. 

A further object of the present invention is to provide a method for calculation of 
crosses between strobe tracks including target type information and the quality of 
crosses, in target type space. 

Still a further object of the present invention is to provide a method for multi-sensor 
multi-target tracking including target type information. 

The object of the present invention is achieved by a process exhibiting the features 
set forth in the claims. The process of the invention uses time recursive Bayesian. 
methods for calculating probabilities for different target types, using different types 
of discrete information, such as ESM data, IRST information and direct 
observations. In order to keep a low complexity in the calculation, approximations 
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are introduced in the likelihood calculations, and an ambiguity restoring procedure is 
introduced for removing inherent abnormalities in Bayesian methods. 

Target type estimates can be used during different stages in a target tracking 
process, e.g. in the strobe tracks and in the target tracks, particularly for improving 
the data to track association, cross calculations, cross quality calculations, track 
quality calculations and for supporting a MHT (Multiple Hypothesis Tracking) 
functionality. When crosses between strobe tracks are calculated, the target type 
information can improve the quality evaluation of the crosses, helping suppressing 
"ghost targets", i.e. false crosses. The inclusion of target type information in the 
track state can also improve the data to track association. When the association of 
data to tracks is ambiguous, several alternative hypotheses may be retained and the 
decision postponed until further data become available, which MHT functionality can 
be improved if also target type information is available. 

Drawings 

Some explanatory embodiments according to the present invention are presented in 
detail in the following, in connection with the associated drawing, in which: 

FIGURE 1 is an example of a possible target type system; 

FIGURE 2 is a block diagram of a preferred embodiment of a target type 

estimation system according to the present invention; 
FIGURE 3 is another example of a possible target type system, including an 

ambiguity possibility; 

FIGURE 4 is a diagram illustrating a simulation of a real target type estimation 

according to the invention; 
FIGURE 5 is another example of a possible target type system, illustrating the 

use of IRST information; 
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FIGURE 6 is a block diagram of a target tracking system according to the state of 
the art; 

FIGURE 7a and b illustrate track association situations, where discrete 

information is used in order to avoid erroneous associations; 
FIGURE 8 illustrates a calculation of a cross from two strobe tracks comprising 

target type estimation information; 
FIGURE 9a and b illustrate cross quality calculations, comprising target type 

estimation information; and 
FIGURE 10a and b illustrate a situation, where target type estimation information is 

used for excluding "ghost targets". 

Detailed description 

A tracking process is a process in which measurements from active and passive 
radar systems or the like are used to estimate mainly kinematic quantities for a flying 
object. Plots, i.e. measurements stating the direction to (azimuthal angle and 
possibly elevation angle) and the distance between a flying object and the sensor, 
as well as strobes, i.e. measurements stating only the direction to a flying object 
(azimuthal angle and possibly elevation angle), are used as input data. 

In a general tracking system, a system track is normally generated in some stage of 
the process. The system track is typically determined by a number of estimated 
kinematic quantities, such as track number, position, velocity, accelerations etc., 
which are generally results of filtering processes. These estimated quantities are 
normally accompanied by uncertainty measures and correlations in the form of 
covariance matrices. The system track may thus be represented as a vector X and 
a covariance matrix P. The system track contains information on where the target-is, 
what the speed is and where it is heading. System track will henceforward be used 
in order to denote such a unique quantity, estimated in a target tracking system. A 
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target track or a track denotes a more general quantity that may exist in different 
parts of the tracking system. 

Note, however, that there is most often no information on what type of target the 
track represents. Target type will in this context mean category or specific type of 
the actual aircraft in question. There may also exist other type of tracks during a 
tracking process, such as strobe tracks. 

For each target type system, i.e. for each situation where a target type estimation is 
supposed to be performed, a predefined number fj of possible target types 7~ are 
available, see e.g. fig. 1: 

T x means that the target is of type 1] , while {7;}^ represents the set of all possible 
target types. 

Measurements are then performed, and different discrete (enumerable) results are 
possible to obtain. Such results - measurement outcomes - are denoted M ' : 

M J means that the measurement outcome M } is detected, while 

represents the set of all measurement outcomes. n m is then total number or 
measurement outcomes in the target type system. In a general situation each target 
is associated with at least one measurement outcome. In the same manner, at least 
one target is associated with each measurement outcome. An example of a possible 
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target type system is shown in fig. 1 . in this example, three target types and four 
measurement outcomes are available. Target type 7;, may give rise to all 
measurement outcomes M x to M 4 , target type T 2 may give measurement outcomes 
M ) , M 2 and M l% while target type T 3 may only give rise to measurement outcomes 
M } and M 2 . It is thus understood that a measurement detecting measurement 
outcome ^indicates that target type 7j is present, while the measurement 
outcome A^may origin from one of the target types 7; or T nt only excluding T z . 
Such information of the relations between target type and measurement outcomes 
have to be added to the target type estimation process based on general 
considerations. The relations are possible to determine by a matrix A™ of 
dimension N T *N M , where "1° denotes a possible relation and "0" an impossible 
one. A matrix corresponding to the example of fig. 1 therefore is, with target types 
running along the rows and measurement outcomes along the columns: 



1111 
1 1 1 0 
110 0 



Target type estimates are added to the system track by including probabilities for 
the different target types in the system track. Unlike a continuous estimated quantity 
and its associated uncertainty, like X and P % the target type is represented by a set 
of probabilities, one for each target type in the limited group of possible target types, 
i.e they may be represented by a vector of discrete probabilities: 



/>(4 #' = !,...,#,. 

The probabilities are estimated based on a series of measurements, giving one of 
the measurement outcomes as a result. The measurements are denoted: 
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■ ■*.). 

where t q denotes the measurement time and where is an index running from 1 
and up indicating the number of the measurement. The total set of measurements, 

<<<)=[i'MU •--(<.)!• 

is the available information for creating the probabilities. 

The use of the measurements can be determined from Baye's rule and the total 
probability theorem as follows, (see e.g. Y. Bar-Shalom and X.-R. Li, "Estimation 
and Tracking: Principles Techniques and Software", Artech House, MA, USA, 1993, 
chapter 1 ) 

for /= I JV r (1). 

That means that the probability for one particular target type at a certain time i q , 

given all measurements, is equal to the probability for the same target type, given all 
previous measurements and given the latest measurement. The probability can then 
be expanded into the normalised product of two factors. The first of these is the 
likelihood of obtaining the last measurement, given the target type and all previous 
measurements. The second factor is the probability for the particular target type at 
the present time but using only the previous measurements. The second factor is 
hence propagated (predicted) to the present time after the last update. This 
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propagation is described below. Initial values are needed in order to start the 
recursion. They are given by the user selected prior target type probabilities 

^oJ,\Z0 0 ))=P,\ for / = 

After a measurement has been used for updating the target type probabilities, the 
target type probabilities need to be propagated to the next measurement time for a 
new update and so forth. The propagation rests upon the following facts. First it is 
noted that target type information is reduced in quality between measurements (as 
is all other probabilistic information). Secondly, it is desirable that the effect of 
measurement errors decay as the time from the error increases. This is achieved by 
letting the target type probabilities decay exponentially towards the prior target type 
probabilities between measurements 

/>(/, , 7;|2(/, rl )) = />° + (/>( Vl , t;|Z( Vl )) - p;)e' ( '<-'<^ 

for / = l,...,A' r , 

where r denotes the decay time. It is obvious to anyone skilled in the art that other 
decay relationships than an exponential can be used as well. 

In order to complete the target type estimator, the likelihood has to be evaluated. It 
is understood that without approximations the computational task rapidly exceeds 
the practical limits, since the complete measurement history is involved in the 
calculation at each step. Approximations are therefore introduced. Examples will be 
described in detail further below. 

An inherent problem in approximate and exact Bayesian methods is the inability to 
accurately handle ambiguities and ignorance (see e.g. S.S. Blackman "Multiple- 
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Target Tracking with Radar'Applications", Artech House, MA, USA 1986, p. 381). 
This will be focused on by an example. 

Example 1 

in fig. 3, a very simple system is shown. Two possible target types T x , 7, and two 
possible measurement outcomes M ]t M 2 are available. The measurement outcome 
M s is only possible for target type 7], while measurement outcome M 2 is possible 
for both target types. The above calculations are first approximated by only retaining 
the latest measurement outcome in the likelihood and by excluding the propagation 
( T - oo ). Then 



K " 3 ' ' .^,.r(^)|7;)^v 1 ,7;jz(v,o) + ^/».^)|7':Hv,,7;|Z(v,)) 



Assume that the system continuously measures M 2 . This means that either of the 
target types may be present, and the probability of the two types should not change 
from the prior target type probabilities, since there is an ambiguity in the target type 
determination. The approximated likelihoods become (see fig. 3 and note that all of 
the measurement history is now neglected) 



p(< q ,M y \T) = L u 



y(/„Hir : )=o 
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II 



P(i,M 2 \T 2 ) = 



Here it is assumed that there are no measurement errors. Note that it is evident from 
fig. 3 that L n +L 2I =K Furthermore, 0<L M <1 and 0<L 2i <] or at least one 
possibility in fig. 3 could be removed. Since M 2 is measured, the target type 
estimation algorithm becomes 

I^P{l q . i ^Z{, q . i )) + \P( t ^J,\Z{i <l .,)) 
Dividing these equations gives 



Since 0< Z, 2I < l it is clear that <9(^) — ^ > 0, which means that 



i.e. a decision is made despite the fact that the situation is ambiguous. An artificial 
instability hence occurs. It can be shown that in cases where the target types have 
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.different numbers of possible measurement outcomes associated with them, this 
instability occurs. 

One of the important features of the present invention is that the artificial instability, 
described in example 1 above, is avoided by modifying the method so that the 
estimated probabilities retain ambiguities in situations where such exist. This 
ambiguity restoring procedure is described in detail further below. A further 
important feature of the present invention is that the computational load is kept low 
since approximations in the likelihood calculations become possible with the 
ambiguity restoring procedure. 

One of the simplest schemes possible for restoring the probabilities in situations 
with ambiguity would be to ignore any measurement outcomes, giving rise to such 
an ambiguity. The probabilities will of course not be changed, but information will be 
lost. In the tracking situation illustrated in fig. 1, the measurement outcome M 4 will 
be used in the target type estimation. All other measurement outcomes M i to M 3 
will be ignored. It is obvious that such an ambiguity restoring procedure will result in 
some improvement but it does not operate efficiently in general. 

In a prefered embodiment of the present invention the ambiguity restoring procedure 
is realised within the scheme of approximation. Subsequently, such a preferable 
embodiment is described in the following. A first approximation is that only a certain 
number of the latest measurements are used in calculating the likelihood, but at 
least one per track. Generally, this is expressed as follows 
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, In the following, the combination of approximations to obtain a low computational 
load and the ambiguity restoring procedure is illustrated for the case where the two 
latest measurements are retained, i.e. A? = 2 . It is a feature of the invention that the 
use of more than the latest measurement allows also good suppression of random 
disturbances. It would be clear for anyone skilled in the art how to generalize the 
procedure to Aq>2 and to apply it for A? = 1 . 

With the two last measurements available, the information thus consists of the 
measurement outcomes 

*(',) = M n 
*(',-. ) = ^,-,- 

Using standard relations (see e.g.Y. Bar-Shalom and X.-R. Li, "Estimation and 
Tracking: Principles Techniques and Software", Artech House, MA, USA, 1993, 
chapter 1) the likelihood can be rewritten as 



' i\t<J,.Z{t q _ x )) />(/,, 7;,z(v,)) 

= P[' q , )l T t , *(/,_, ). Z(/,_, ))P[t, . z(/,_, )| T, , Z( V : ) ( , 

/ = 1 N T . 



Then consider the quotient in the right hand side of the last equality. When the 
effect of propagation is small and when more and more data are collected, the 
numerator and the denominator will be using almost the same information. In the 
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flmit where q-^co this will hold true exactly as formulated in the following 
approximation 

This reduces the likelihood to the product of two factors 
i = \,...,N T . 

The final approximation is now to truncate the measurement history according to 
Aq = 2. This results in the following likelihood 

^/,,r(/,)|7;,Z(^)).p(/,,AY,|7;.W f . 1 )^/,.^ l |7;) 1 i = \,...,N T 

The above result can be represented as a look-up table, from which the requested 
probability values are picked when needed for the probability calculations. If only 
the two latest measurement outcomes are used, the look-up table will consist of a 3- 
dimensional data structure of dimensions N T x N M x N M . In a real system, the 
number of possible target types could typically be of the order of 25-500 and the 
measurement outcomes could be several thousands, why the data structure is 
usually large. It is obvious to anyone skilled in the art that it is also difficult and 
cumbersome to find all conditional probability values. In the preferred embodiment 
of the invention the look-up table is constructed in such a way that the mpst 
important effects are captured by a small number of parameters. At the same time 
the ambiguity restoring procedure is introduced. One such possible scheme 
according to a prefered embodiment is presented next. 



WO 98/58274 



PCT/SE98/00956 



15 

Above, it was shown that it is advantageous to introduce an ambiguity restoring 
procedure. The origin of the problem was shown to be the fact that different target 
types can have differerent numbers of measurement outcomes associated with 
them. The first major step is then expressed in the following statement: 

Define one of the target types to be the designing target type, denoted by T D . 
T n is used for computation of those conditional probabilities that are related 
to ambiguities (described below). The obtained values are then applied to all 
the other target types using a special procedure. The other conditional 

probabilities, that are not related to ambiguities, do not need consideration of 
T D in the design. 

With reference to example 1 above, this means that if T D is chosen equal to T 2 , the 
modified result would become 

p[> q ,M i \T)=L u 

The effect is thus that p[t q ,M 2 \T^ is modified from L,,to 1. As a consequence, 

following the reasoning of the above example, no decision will be made in 
ambiguous situations. 
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In the following, the above described procedure involving a designing target type is 
applied to the case where the two latest measurement outcomes are retained. The 

design of the first factor of the likelihood, P^^M^T^M^) , is first handled. The 

following parameters are used 

P c : the probability of correct detection of measurement outcomes. 

P E : the probability of erroneous detection of measurement outcomes. 

P f : the probability of detection of measurement outcomes that are consistent 
(the measurement outcome is possible for the target type) with the 
possible measurement outcomes, for the designing target type. 

The parameter P E is related to P c as 



Jy M 1 



since there is only one way to do a correct detection in the case without ambiguities. 
The parameter P l follows from the probability equation 



{N M 'N M (d))F E ^N M (J)P I =\ 



since there are N M {d) possible consistent measurement outcomes for the 
designing target types and since the rest are erroneous. Thus 



N M -\-(N M -N M (d))(\-P c ) 
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All parameters can therefore be computed from knowledge of the probability of 
correct detection for the sensor. This single detection probability therefore 
constitutes the only input parameter, except the original definition of the target type 
system, which is necessary to define the probabilities in the look-up table. It is 
important to note that these three parameters, originating from the same detection 
probability, are used for all target types in the scenario. One set of parameters thus 
results for each sensor 

The above three parameters are then used in order to compute the look-up table 
corresponding to p[i q ,M q \T t ,M q ^. This is done by consideration of the following 
four cases 

1) M q -M q ^ and M q _ x consistent with 7;. 

2) M q *M q _ x and A^_, consistent with T { . 

3) M q ^ inconsistent with T it ^consistent with T { . 
A) M q _, inconsistent with 7;, M q inconsistent with T f . 

These cases covers exactly all possible situations when the two last measurement 
outcomes are retained in the likelihood. It is assumed that all measurement errors 
are random. 

1 ) In this case the two consecutive measurement outcomes are equal and 
consistent. Because of the high reliability (same measurement outcome and 
consistent measurement outcomes) this case is assigned the probability P c . 

2) In this case the previous measurement outcome was consistent with T t while the 
present one is different from the previous. The difference indicates a 
measurement error since the previous measurement outcome was consistent. 
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Measurement errors can also give rise to consistent or inconsistent measurement 
outcomes. Hence this case is assigned the probability P E . 

3) In this case the previous measurement outcome was inconsistent with 7,', which 
indicates an error. Then there is no information in the previous measurement 
outcome when the probability of the present measurement outcome shall be 
determined . In case 3) the present measurement outcome is consistent with T t 
which has the probability P } associated with it (neglecting the measurement 
history). 

4) In this case the previous measurement outcome was inconsistent with T t t which 
indicates an error. Then there is no information in the previous measurement 
outcome when the probability of the present measurement outcome shall be 
determined. In case 4) the present measurement outcome is inconsistent with T t 
which again indicates a measurement error. Hence the probability P E is 
associated with case 4). 

By using the target type system definition matrix A™ , the probabilities can be 
expressed as: 

+ P £ (A-. i+ (l-A™).(l-A^)) (4), 

where / = l,...,tf r , N M% A^_ } = 1,...,A^ and 

S„ Ai is the Kronecker delta function. 

For the second factor of the likelihood, A/,.,17;), it is sufficient to distinguish 
between measurement outcomes consistent with T that are given the probability P } , 
and measurement outcomes that are inconsistent with T t that are given the 
probability P E . 
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In terms of indices, the relation is: 

P('<,M^T) = P E +(P,-P g ).A™^ (5), 
where /=1,. A^_, = \,...,N M . 

This completes the likelihood calculation. 

The simplest example of discrete measurements are direct observations. In such a 
case, there is a one-to-one agreement between target type and measurement 
outcome, why the target type system in such a case is defined by a square matrix 
with T in the diagonal. 

Example 2 

A more relevant example of the above method is given by the situation defined by 
fig. 1. The measurements are in this example ESM measurements and accordingly 
the measurement outcomes are emitter modes. The look-up table is thus a 3- 
dimensional data structure with dimension 4x3x4 and a 4 x 3 matrix. In order to 
visualise the data structure, each target type is illustrated separately, i.e. 3 separate 
matrixes 4x4 together with the 4 x 3 matrix. Target type 2 is selected as the 
designing target type. According to (2) and (3) the following parameters are 
obtained: 

P c : user choice. 
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The data structures become 



p(> q M,\T„M q _,) = 



Pc Pb Pb Pb 

Pe Pc Pe 

Pb Pb Pc Pb 

Pb Pb P E Pc 



p p 



E 

Pc 



p, 
p, 



p{l q ,M q \T„M q _) = 



' E 

Pb 



Pb 
Pc 



E 

Pb 



P, 
P, 



P, 



Pr Pr P, 



E 



>, P; P,' 

P; P> P, 

P, P> Pe 

f> P E Pe. 



according to (4) and (5). In the three first equations above, M^runs over rows and 
M q _ x runs over columns. In the last equation above runs over rows while T t 

runs over columns. Starting with p{t q . M q \T^M q _^< a previous measurement 

outcome which is consistent with the considered target type T x and which is identical 
to the last measurement outcome will have probability P c . Since all measurement 
outcomes are consistent with the target type 7J, the matrix has P c along the 
diagonal. The elements beside the diagonal correspond to the cases where both 
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measurement outcomes are consistent with the target type 7; (which in this example 
is true for any measurement outcome), but not equal. The new measurement is then 
considered as an error measurement with probability P E . The whole matrix is thus 
fully determined by the sole factor P c . 



The matrix P\t^M q \T^M n _ x ) is built up slightly differently. Again, previous 
measurement outcomes which are consistent with T 2 and identical with the new 
measurement will be set to P c . Since T 2 is consistent with M ] , M 2 , and M 3 the 
first three diagonal elements are equal to P c . In the same way as above off- 
diagonal elements except in column 4 are set to P E . The elements in column 4 are 
connected to situations where the previous measurement was not consistent with 
the actual target type. However, if both measurements are inconsistent, in this 
example measurement outcome M A , the result is considered as an error and will be 
assigned to the same probability as the other error measurements. The remaining 
elements correspond to cases where the previous measurement was inconsistent, 
but where the new one is consistent. These measurement outcomes are assigned 
the probability P } . 

The matrix p(( q M q \T 3 , M ^) is built similarly. Rows 1 and 2 are filled with P c al the 
diagonal and /^outside, since T 3 is consistent with M x and M 2 . Columns 3 and 4 
corresponds to situations with inconsistent previous measurement outcomes. 
Finally, the 4 x 3 matrix P{j q .M q _ y \T^ contains P } whenever the measurement 
outcome is consistent with the target type and P £ otherwise. 

In this way the whole data structure, i.e. the look-up tables, are defined using only 
the target type - measurement outcome definition A m and one parameter (e.g. P c ) 
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as shown by (2)-(5). By using the concept of a designing target type t all matrices are 
using the same set of probability values and thus any ambiguity is restored. 

Simulation 1 

The above example 2 is in the following used together with a series of 
measurements. The a priori probabilities for the three target types are equal, i.e. 
1/3. The set of true emitter mode transmissions consists of a 100 totally, one every 
tenth second; 25 of measurement outcome M ] followed by 25 of measurement 
outcome M 2t 25 of measurement outcome A/ 3 and finally 25 of measurement 
outcome M 4 . The true transmissions are corrupted with a statistical error as high as 
30%. Someone skilled in the art immediately concludes that 7; is the correct target 
type. P c , i.e. the probability for a correct measurement, is set to 0.70 which is in 
agreement with the simulated measurements. The obtained target type estimates 
are illustrated in fig. 4. It can be seen in that figure that the scheme operates well 
despite the large amount of measurement errors. The three curves P 1( P 2 and P 3 
shows the estimated probability for the target types T K , T 2 and 7^ , respectively. After 
500 seconds (50 measurements), information making it possible to exclude T 3 
becomes available and T y is excluded after some time. The exclusion is soft 
because of the measurement errors. After 750 seconds (75 measurements) 
information becomes available, that allows also an exclusion also of T : , which is 
automatically performed. This is a typical illustration of how the scheme operates by 
excluding target types that are not consistent with the obtained measurements given 
the prior information of possibilities. 

Another example of measurements which are possible to use within the present 
invention are e.g. IRST measurements. The IRST measurements can give 
information about the existence of afterburner on an (air)target by measuring a 
change in intensity, the number of engines on a target, the size of a target (if range 
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is available), missile firings from a target and kill assessments. From such" 
measurements a set of measurement outcomes is definable. 

Example 3 

In fig. 5 an example of a target type system using IRST measurements is shown. 
Four target types are possible, 7J -fighter, T : -bomber, r 3 -fighter bomber and T r 
reconnaissance aircraft. The fighter can carry 6 missiles, the bomber is not 
equipped with missiles, the fighter bomber can carry 4 missiles and the 
reconnaissance aircraft can carry 2 missiles. Accordingly, 7 measurement outcomes 
are possible; M cr 0 missiles fired from the track, M x - 1 missile fired from the track, 
M 2 - 2 missiles fired from the track, M y - 3 missiles fired from the track, M. - A 
missiles fired from the track, M s - 5 missiles fired from the track, hA t - 6 missiles fired 
from the track. The relation between target types and measurement outcomes is in 
this case given by the matrix: 

"] 1 I 1 I 1 f 

10 0 0 0 0 0 

A = 

1 '1 1 1 ] 0 0 

1 1 10 0 0 0 

The approximation scheme described above can then be used. 

in a preferred embodiment of the present invention, discrete measurements of 
different types are used together. If there are measurements available from direct 
observations as well as from ESM and IRST measurements, they are possible to 
integrate in the calculation of one set of probabilities for the different target types, for 
each track. In this way the different data types each give their contribution to the 
exclusion of target types that are inconsistent with measurements and with prior 
knowledge of measurement outcomes and target types, as defined by matrices A 7 *' . 
This is one of the very important advantages with the present invention. The set of 
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measurements then consists of measurements of different types. The invention is 
not limited to direct observations, ESM information and IRST measurements. It is 
obvious to anyone skilled in the art how to integrate also e.g. COMINT 
(Communication Intelligence) information which is closely related to ESM 
information. 

To describe how this is implemented it is most instructive to return to the definition 
of measurements and then to consider the update equation. In this description, only 
two types of different measurement outcomes - type 1 and type 2 - will be discussed. 
The inclusion of other data sources by definition of e.g. measurement outcomes - 
type 3 follows a similar pattern and it is obvious for anyone skilled in the art how to 
proceed. The measurement outcomes - type 1 and measurement outcomes - type 2 
may e.g. correspond to ESM measurements and IRST measurements, respectively. 

The following quantities are needed: 

{7;}^ : the set of possible target types. 





z m(*i) : the present measurement - type 1, measured at time 



2 



('«) : the present measurement - type 2, measured at time t m . 



z{t q ) ; the present measurement - any type, measured at time t q . 
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Z(t q )= : the set of all available measurements 



Note that the subscripts , and m do not take all values since they form a partition of 



The purpose is to be able to use the same methods as defined above. In the 
previously disclosed process, a key detail is that the interdependence of different 
data types does not pose a problem, since there is only one type of data. The 
following simplifications are therefore used in the calculation of the likelihoods, 
when data from different measurement types are present: 



P(t< yMW* ZO, )) - K'« ■ r « (/ < }l * • Z(/ < >) = • z " (/ ' }l T ' ■ Z( '' } ) • = »• * r 



This means that each data type is handled separately in the calculation of the 
likelihoods. The likelihoods are, however, used for updating of one single set of 
target type probabilities. The previously described methods can with this 
simplification now be used for simplified computation of the likelhoods of the 



the subscript 9 that does take all values. 



The starting point is the previously discussed update equation, c.f (1) 




1=1 
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different data types. A very important feature with the method of the present 
invention is that this fusion of different sets of discrete data is easily carried out. 

The probability calculations according to a prefered embodiment of the present 
invention follow in general terms the procedure illustrated in fig. 2. Here, the flow of 
the different measurements are shown together with the dependence of the user 
selected probability of correct detection and the target type system definition. A 
measurement is performed and the result is sent to box 10. In this box , the 
measurement will be delayed one step in order to become the second latest 
measurement M q ^ when passed further to box 1 1 . The latest measurement will also 

pass directly to box 11, then being the latest measurement M q . The probability for 

correct detection P c and the definition matrix A™ of the target type system is used 
in box 12 to compute the likelihood with an ambiguity restoring procedure to form a 
likelihood look-up table. The likelihood is computed in box 11 using the look-up 
table and the two last measurements M q , M^ x and is used to update the target 

type probabilities in box 13. The updated type probabilities are in box 14 propagated 
to the time for the next measurement / t which originates from the measurements, 

and propagated type probabilities are delivered back to be used in the next update 
of target type probabilities in box 13. 

As mentioned in the introduction of the detailed description, type probabilities are 
possible to use in different stages of a target tracking procedure. In the following, a 
preferred multi-sensor multi-target tracking system, which in part is disclosed in the 
Swedish patent application No. 9700709-0, is summarized as a background for other 
possible uses for the present invention, but first a few definitions are needed. 

A plot is defined as a direction and a range, i.e. a position in space. Range 
information is available if the sensor is of an active type. A measurement from a 
passive sensor has no range information, why the results are purely bearing 
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information, and in this case the measurement is called a strobe. If several strobes 
from the same sensor are assumed to origin from the same target, a strobe track 
may be initiated, i.e. a track in the angular coordinates with reference to a certain 
sensor. Where two strobes or strobe tracks crosses, strobe crosses and strobe track 
crosses, respectively, may be calculated. These strobe crosses or strobe track 
crosses correspond in some sense to plots, since they define a position in real 
space. A quantity, cross quality, can be defined for every cross, which indicates the 
magnitude of agreement between the involved strobes or strobe tracks. A cross is 
said to be of 2 nd order if two strobe tracks are used to form it, of 3 rd order if three 
strobe tracks are used, and so on. 

The preferred multi-target tracking system is best described with reference to fig. 6. 
Note that strobes and tracks will generally carry type related information when the 
tracker is upgraded with the invention described herein. 

Plots 20 and strobes 21 are first preprocessed in an input processing unit 22 and 
collected in time batches in a scheduler 23. The plots 20 and the strobes 21 are 
then sent to a tracking kernel 24 for association to existing system tracks 32. 
Associated plots 20 and strobes 21 are used for updating of these system tracks 32. 
Plots 20 that are not associated to system tracks 32 are sent to a plot handler 25 for 
initiation of new system tracks 26. Strobes 21 that are not associated to system 
tracks 32 are sent to a strobe handler 27 for initiation of strobe tracks, calculation of 
strobe track crosses 29 and finally initiation of new system tracks 26. Bayesian 
algorithms are used for triangulation (creation of strobe track crosses 29) and 
deghosting. A bias calculator 30 continuously computes bias compensation 
parameters 28 from measurements that are firmly and uniquely associated to high 
quality tracks. Such bias compensation parameters 28 are stored in a sensor 
database 31 and used during the tracking process. .System tracks and strobe track 
crosses 29 are sent for display. 
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The MST (Multi-Sensor Tracker) makes full use of sensor oriented coordinate 
systems for maximum flexibility. This makes it possible to handle an arbitrary mix of 
active and passive measurements at all times. 

All system tracks are represented and propagated in a cartesian, 3-dimensional, 
earth tangential coordinate system, centered somewhere in the surveillance area. 
The system tracks are transformed to map coordinates for display. The association 
of data to system tracks is performed in the measurement space to allow a handling 
of passive measurements without range information in a coherent manner. The 
system tracks are thus transformed to (possibly a lower dimensional) measurement 
space where a conventional multiple hypothesis maximum likelihood association is 
performed, see e.g. S. S. Blackman "Multiple-Target Tracking with Radar 
Applications" , Artech House, MA, USA, 1986, pp 249-280. To improve the 
performance, scan to scan memory is introduced in the association process by the 
use of multiple hypothesis tracking (MHT) (S. S. Blackman, Multiple-Target Tracking 
with Radar Applications", Artech House, MA, USA, 1986, pp 289-300). The MST 
described here implements a variant of MHT that is sometimes denoted track 
oriented MHT. All associations with a sufficiently high quality, together with a 
coasting alternative, are then retained. The resulting hypotheses are used in order 
to form alternative tracks, for each system track. The probability that each 
alternative is correct, conditioned on the fact that there is only one target, is updated 
with Bayesian techniques. Hypotheses are removed when they fall below a 
threshold. 

The system tracks are also updated in the measurement space of each sensor, 
using extended Kalman filter (EKF) techniques, see e.g. Y. Bar-Shalom and X. -R. 
Li, "Estimation and Tracking: Principles, Techniques and Software", Artech House, 
MA, USA, 1993, pp 382-399. 
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The initiation of tracks is treated separately for active (plots) and passive (strobes) 
measurements. The active initiation relies on multiple hypothesis techniques, from 
scan to scan. The automatic passive initiation is much more elaborated. Briefly, it 
relies on triangulation of strobe tracks to form crosses. The quality of each cross is 
then evaluated with very advanced Bayesian techniques. Adaptive clutter density 
maps are automatically updated in the system to support the track initiation process. 

The obvious objective for the target type identification method achieved by the 
present invention is thus to be used to identify a target type for each track (strobe 
track, target track or system track) in a tracking system as the one described above. 
The additional information that is created by estimating target types can also be 
used in order to improve existing functionality in a kinematic tracking system as the 
one described above. 

To describe this, each (kinematic) track (strobe track or target track or system track) 
of the system is first extended with a discrete part, said part consisting of the 
estimated target type probabilities as described above. This discrete track is 
updated and propagated according to the descriptions above. 

A necessary prerequisite for using the target type estimates in other parts of the 
tracking system is the computation of quantity denoted the discrete measurement 
probability function. This quantity is the probability of the measurement (discrete set 
of measurement outcomes) given the estimated target type probabilities of a 
particular track. The computation is technically involved and the details can be 
found in appendix A. The case where the two latest measurement outcomes are 
retained is discussed here. The result is 




1=1 
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Here, the index p denotes measurement number p in a set of measurements. The 
index n means that the quantity belongs/refers to (discrete) track number n . The 
quantity / D (/,,/>,«) is therefore the probability of the measurement p given all 

information condensed in track n. This quantity thus measures how well the 
particular measurement correlates to the particular track. The result of the 
computations is thus a sum over quantities consisting of three factors. It can be 
seen that the two first factors correspond to the likelihood discussed previously. 
However, as previous measurement the previous measurement that was used for 
updating of the particular track shall be used. The last factor consist of the 
estimated target type probabilities of track 

Association ■ 

The discrete probability functions achieved by the above described method are 
possible to integrate with continuous probability functions based on the kinematic 
parameters, e.g. in the target tracks achieved by the above described process. To 
describe this process it is suitable to first discuss the basics of data to track 
association in kinematic multiple target tracking and to illustrate the use of discrete 
information in the association process. It is obvious to anyone skilled in the art that 
the kinematic association process described here is one of many possibilities and 
that the present invention regarding the integration of target type information may be 
applied to other kinematic association processes than the one described herein. 

When the system receives measurements from the sensors it collects them in time 
sliced batches. As soon as a batch is filled with measurement reports it is sent to the 
tracking system. Before the measurements can be used to update the present tracks 
in the tracking system it must be determined which measurement (if any) that 
corresponds to a specific track. This procedure is usually denoted association or 
correlation. In one possible implementation the procedure can be roughly divided in 
three main steps. 
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Step 1: The first step is to exclude tracks that are out of range of the measurements 
of one particular sensor in the present time slice. This means that min/max limits of 
range, azimuth and elevation are computed and that those tracks that do not fall 
within these limits are excluded from the association procedure. In the case where 
discrete information is available one can also use known limitations of for example 
ESM sensors to exclude tracks that cannot possibly have generated the 
measurement. For example, band limitations of sensors may exclude emitter mode 
measurements in certain frequency bands, which may in turn exclude tracks 
associated with particular target types. The tracks that fulfil the limitations of the 
sensor are the only tracks that are considered in the subsequent association 
procedure. Step 1 ends by propagation of each remaining track to the mean time of 
the remaining measurements. (This should be close to the mean time of the time 
slice provided that there is a significant number of measurements.) 

Step 2: The second step in the measurement to track association comprises forming 
measurement and track clusters. A cluster is a set of measurements and tracks such 
that 1 ) The measurements are exactly those that are included in the union of the 
track search regions and 2)The clusters cannot be divided into smaller clusters. 

It is important to note that the search regions are defined in the measurement 
space. This space can include kinematic dimensions as well as dimensions 
corresponding to measured discrete information. For example, the measurements 
may consist of azimuth and emitter mode from an ESM sensor in which case the 
measurement space has two dimensions (continuous azimuth and discrete type). 
Since the dimension of the track may be higher than that of the measurement space 
the track is transformed to the measurement space. 

The clustering procedure starts with a measurement (a point in the measurement 
space). Tracks that are within a search region centered around the measurement 
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- are. then included in the cluster. For each track in the cluster, measurements that 
are within a search region centered around the track are then included in the 
cluster. The procedure continues until it terminates or until a maximum number of 
measurements and tracks are obtained for the particular cluster. In order to make 
this procedure successful, the search regions need to be increased from pure 
measurement and track uncertainties to account also for timing errors. These errors 
occur because the mean time of the measurements must be used in the clustering. 

Step 3: The third and last step is a multi hypotheses test procedure based on 
statistical maximum likelihood techniques. The purpose is to find an optimal 
association of measurements to tracks within each cluster. The procedure is 
denoted the LPQ (Logarithmic Probability Quotient) and it is described in some 
detail below. 

To explain the LPQ method it is assumed that there is precisely one measurement 
considered and the task is to determine whether this measurement shall be 
associated with the track in question. When doing so the following two constraints 
must be fulfilled. 

1 ) At most one measurement can be associated with each track. 

2) Each track generates at most one measurement. 

The following hypotheses are now relevant . 

H x ) The measurement is true, originates from a true target and no false 

measurement is obtained. 
H 7 ) The measurement is false and originates from a false target. 
H y ) The measurement is false, originates from a true target and this target is not 

detected. 
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Jt is also assumed that the probabilty of detection of the sensor is P Dt that the 
probability that the measurement is spurious (false) is P s . and that the probability 
that the track corresponds to a true target (probability of true target) is P(7T). 
Furthermore, the probability density, function (pdf) of the kinematic part of the 
measurement conditioned on the predicted track is denoted a(^,/v?). This pdf is 

usually Gaussian and it may be evaluated in the track coordinate system (3-D 
measurements) or in the measurement space. The probability density function of the 
discrete part of the measurement conditioned on the target type probabilities of the 

track and previous measurements, is denoted f D [i q ,p,n). 

The next step is then to consider a cluster of measurements and tracks. The task is 
to determine which measurement (if any) that belongs to a particular track. This 
requires a measure that determines the quality of a tentative association of track n 
to measurement p (cf. the discussion above on the discrete probability functions). 
Towards that end, consider the hypotheses above and assume that all phenomena 
above can be treated as statistically independent. Then the likelihoods of the 
hypotheses above can be expressed as follows 

pit^H^p^ccP^TDf^i^p^foit^p^) 

since normally P s *0. Here, the kinematic part of measurement p is denoted 
X n {i \i x ) denotes the propagated kinematic state of the track and 
P" (t q \t denotes the corresponding covariance matrix (of track n). Note that this 
probability density function is defined in the measurement space and that it is 
normally Gaussian. Considerations similar to those above give 
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p{t q M„p,n)ozP s (\-P(TT)) 



p(t q ,H 3 ,p,n)ccP s P(7T)(\-P D ) 



As a measure of whether it is good or bad to associate measurement p with track 
n it is now very natural to use 



The reason is that under // 2 and // 3 there is no reason to do an association, while 
this should be considered an alternative under //,. The logarithmic probability 
quotient ( Ipq ) of these events is then a natural score function. Note that the score 
function will be positive if is more likely than the other hypotheses when 
combining track n with measurement p. Some algebra results in 



Having computed the score function for a general track to measurement association 
the best combination among all possible such associations need to be determined 
since this solves the problem. Towards this end the following matrix is introduced 
where the rows are indexed by measurements and the columns by tracks. 





(P D P{TT)f K {i^p,n)f 0 {i q ,pji) 



P S (\-P D P(TT)) J 



LPO = 



lpq{t qP ,P,\) lpqQ qP ,P y N) 
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The optimization problem now consists of choosing at most one element of each row 
and at most one element of each column so that the sum of the chosen elements is 
maximized. Methods that solve this problem is readily avalable from e.g. S. S. 
Blackman, "Multiple-Target Tracking with Radar Applications", Artech House, MA, 
USA, 1986, pp 397-401. 

The use of the discrete measurements in the last step of the association process 
can thus be summarized as follows: 

Prerequisites: 

1) Each track is extended with a discrete state consisting of target type probabilities. 
Each track is updated and propagated according to the methods outlined above. 

Data to track association: 

1) Compute the discrete probability function for any measurement to track 
association. 

2) Replace the kinematic only probability density function f k (t qy p,n) with the 
product f K (t qi p J n)f D (t qJ p 1 n) in the computation of lpq{t q% pji). 

3) Run all existing functionality for kinematic data to track association as usual. 

The association process has above been described for tracks. Anyone skilled in the 
art will understand that this process can be extended to different kinds of tracks, 
than the ones discussed herein. 

Example 4 

To illustrate the improvements that can be expected, fig. 7a and 7b show a. scenario 
with two targets T1 and T2 : where the filled shapes symbolises the different target 
types, and two measurements (strobes) S1 and S2 that are to be associated. The 



WO 98/58274, 



PCT/SE98/00956 



36 

line extending from the target track indicates the velocity by the direction of this line. 
In fig. 7a only kinematic information is used, i.e. no significance is paid to the target 
type information, while target type information is used in fig. 7b. The angular 
accuracy's of the strobes are assumed to be significantly larger than the angular 
separation between the strobes. It can be seen that the association in fig. 7a fails 
and associates strobe S1 to target track T1 and strobe S2 to target track T2. 
However, the additional information from the target type estimation, symbolised by 
the empty shapes at the end of the strobes, make a correct association possible in 
fig. 7b, i.e associating strobe S1 to target track T2 and strobe S2 to target track T1. 

The above four step procedure (prerequisites and data to track association) is in 
fact of general validity. An important advantage of the present invention is that it 
uses this procedure in a multitude of existing functions in the kinematic target 
tracking system. The use and benefits of using the procedure in data to track 
association (to system tracks) have been illustrated above. Next, the use of the 
procedure in other parts of the tracking system is discussed. 

Strobe tracking 

In an MST system as the one described above a method for passive initiation of 
multiple tracks is available. The first step is there that strobe tracking is performed 
for each sensor. The process of strobe tracking comprises a full multiple target 
tracking process in the angular domain, that is performed for each sensor. In that 
process strobe tracks (essentially Kalman filter states in the angular domain) are 
maintained and incoming strobes are associated to the strobe tracks. In this process 
discrete target type related data can be used exactly as in the association of strobes 
to system or target tracks, using the four step procedure above. The only difference 
is that the kinematic tracks now reside in the angular domain. 

Computation of track quality 
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.In order to compute the logarithmic probability quotients it is necessary to have an 
estimate oiP(TT). It is the purpose of this section to show how this value is 
computed when discrete information is available. The probability of true target is the 
propagated value of the quantity when used for association purposes, see below. 

The scanning mono sensor case is considered here. In this case it is known whether 
a measurement miss has occurred since it is possible topredict the next time of the 
measurement from a particular track. If no measurement is associated with the track 
for a specific batch and if an association should have taken place during the time 
span of that batch, then a measurement miss has occurred. 

First consider the case when a measurement has been obtained. Now, by definition 



/>(77,/,) = P(TT\Z(t,)) = P(77l2(^),Z(v,)) = 



/>(2(/,)17T,Z(y,))P(771Z(^,,)) 
P(z(t q )\Z(t q .,)) 



P(2(t q )\TT,Z(i^))P(TT,l q . t ) 



P{z{( q )\TT,Z(t q .,))P(TT,t q .,) + P(r(r,)h7T,Z(V l ))(l - W,-)) ' 



Assume that 



P(z(t q )\TT,Z{t^)) * P(z(t q )\TT) and P(z(i q )\-,TT,Z(i q . 1 )) * P(z(t q )\-TT). 



Then it follows immediately that 



P(TT,i q ) = 



PW<)\TT)P{TT, /,.,) + P(z(i q )\-,7T){] - P(TT ,/,.,)) 



p(z{t,)\Tr)PCn\t q _>) 



The likelihoods can then be evaluated. The result is 
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/>(r(/, )\TT) = P D f K (i k , p, ,i)f D (/, ,p, n) + P s ( 1 - P D ) 

P{z(l q )\-,TT)=P S 
This gives the result 

1 

When no measurement has been received when it should the result isinstead, using 
similar approximations as above 



P(-*{t q )\TT,Z{t q _,))P(TTj q „,) 



p s (\ - P D )P(TT,t,) + p s (\ - P{TTj q _,)) : p s \~P(rrj^) 

' \~P D P(TTj^) 



The integration of discrete data is again reduced to a simple multiplication of the 
discrete probability function with the kinematic counterpart. Although the above 
description was for the monosensor case, ir is obvious for annyone skilled in the art 
how to extend it to the multi-sensor case. 

Multiple hypotheses tracking 

As described above, MHT introduces delayed decisions in the association process 
by retaining all ambiguous associations until they are resolved. This is done by 
starting one internal track for each association to the system track. The quality of 
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.these alternative tracks are then continuously evaluated by a track credibility 
measure, as the one described above. Discrete data can be included in this process 
by the four step procedure described above. Alternative tracks with a poor goodness 
of fit in the target type domain will then be more efficiently pruned. 

Calculation of strobe track crosses 

In the above mentioned MST system a procedure for computation of the kinematic 
crosses between strobe tracks is used. In that procedure, strobe tracks (containing 
bearing(s) and bearing rate(s)) are combined to produce a position and a velocity in 
Cartesian coordinates. First a rough gating excludes infeasible combinations of 
strobe tracks. Then all crosses between strobe tracks are computed statistically. 
First, all crosses containing two strobe tracks are computed. Then all crosses with 
three participating strobe tracks are computed and so on up to an order equal to the 
number of sensors that are connected to the data processing center. 

Now strobe tracks can be extended to also contain target type information. This has 
been described in detail above. It is then natural to use the target type information 
also in the computation of strobe track crosses. This type of information in a strobe 
track cross is for example helpful when a system track is initiated from the cross. 
The question is then how to compute the resulting target type state, given the target 
type states of the strobe tracks. This corresponds, in some sense, to the weighted 
averages that occur in the kinematic computation of strobe track crosses. 

Denote the strobe tracktypes by S r j = 1 N s .The sought quantities in the case of 

a strobe track cross of order N are then 



WO 98/58274 % 



PCT/SE98/00956 



40 

.Here the indices j\ to j N is a subset of the available strobe tracks, that are all 
produced by different sensors, T* is here the target type / of the strobe track cross 
X . Baye's rule and the total probability theorem give 



Km *- 

?1P(S*.^SJT I *)P{T?) 



The prior target type probabilities P(T t x ) are obtained from prior information. One 
solution is to set them all equal so that their sum equals 1. The other quantities can 
be evaluated as follows 



, x P[s t . S } , T t x ) Ms , , T t x , S f J* S\ , T* ) 

1 = 1,...,^. 



Assume that the strobe tracks are statistically independent. The following 
approximation is then motivated 

When this is inserted above the final result becomes 
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The quantities P(s jt \T t x )> k ~ 1 N , are the target type probabilities of the strobe 

track, given the information that the target type equals l] x . This is nothing else than 
the target type probability for target type 7/ v , corresponding to the strobe track s . 

This quantity is exactly what is produced by the target type identification scheme, 
when applied to the strobe tracker. The result of the computations is thus that the 
strobe track cross is equipped with a full target type state. 

The final result has a structure where the estimated probabilities of target type 7/ v 
for each participating strobe track are multiplied. This means that if the probability is 
low in one single strobe track it will also be low in the resulting strobe track cross. 
Consequently, only target types that have not been softly excluded in any strobe 
track will have a high probability in the resulting strobe track cross. This effect is 
illustrated in fig. 8. 

In fig. 8, a calculation of a 2 nd order cross is illustrated. The target type probabilities 
are high for target types 2 and 3 for strobe track 1 , and lower for target types 1 and 
4. For strobe track 2, the target type probabilities are high for target type 1 and 2 
and low for types 3 and 4. When the 2 nd order strobe track cross is calculated from 
these strobe tracks, the individual target type probabilities are multiplied. Thus, the 
2 nc order strobe track cross will acheive a high target type probability for target type 
2, since both factors are high. The probabilities for target types 1 and 3 will be 
somewhat lower, since they are formed from one high and one low value, and finally 
target type 4 will be very small, since it is the product of two small factors. In this 
way, it is clear that target type 2 is the main candidate for the strobe track cross. 

Calculation of the quality of strobe track crosses 

The quality of the strobe track cross is a measure of the probability that the cross 
corresponds to a true target, conditioned on the participating strobe tracks in the 
cross. This quantity; among other things, measures how well the participating strobe 
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.tracks fit together. A calculation of the kinematic quality of a strobe track cross of 
arbitrary order is included in the above described MST system. When target type 
information is available, this quality measure can be further refined. Exactly as 
above this is obtained by treating the target type information separately from the 
kinematic information and then using the two resulting measures together to form 
the complete quality measure of the strobe track cross. 

The evaluation of the discrete quality of a cross in discrete space corresponds to the 
evaluation of 

s Ja ) 

which is the probability of a true cross (TX) conditioned on all participating strobe 
tracks used in the computation of said cross ( X ). This quantity can be computed 
recursively over the participating strobe tracks which is shown in the following. First 
apply Baye's rule 

, , piSATXJ, .V, )p(tX\S i S,) 

^ Jl \TX,S h S^TXjS,, S Jt ) 

~ P (SJTX,S, S lK )p(TX\S J ,,...,S jK ) + p(sj-,TX,S h S jK )p(-,TX\S h S Jk ) 

= ^^^^•■••^H^ S /m ) 

" P(s h \TX,S h ,...,S Jm )p(mSj,....^ S,J 

This relation shows how to incorporate the information of one more strobe track into 
a probability of a true cross computed from a number of strobe tracks. To perform 
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the integration the following two quantities need to be computed and this is done 
using approximate methods; 

p(S h HTX,S J: ,...,S jK \ 
The first quantity is approximated as follows 

The dependence on other strobe tracks than the one under consideration ( S h ) is 

thus dropped. To proceed, note that conditioned on the fact that the cross is true, 
the type is one of 7*, / = l,...,N T (which type it is is not known). Hence 

£p(S h \T, x )p{T*) Nr 
p(S h \Tx)=p(s h \T^T^...^) = ^ rr = I4WTM 

(7) 



using Result 2 of Appendix A and the fact that the target types form a mutually 
exclusive and exhaustive set of events. The end result is therefore 

p(s i }TX,S h ,...,S i ^Y J P(s ji \T, x )p{T, x ) (8). 



The computation of the second quantity starts with a similar approximation, i.e 
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Now, the strobe track has been generated in the system from several strobes so it is 
very likely that there is a corresponding target. The probability of the strobe track, 
conditioned on the fact that the cross under evaluation is not true, hence means that 
some of the other crosses on the considered strobe track is the correct one, or that 
the target is seen by one sensor only. Therefore 

^ h7*) = p[S h \TX\ v TX' 2 v...vTX K , v Ts) 

where TX\. , = l A/ denote cross k of all crosses that have the strobe track 

under consideration as a participating strobe track, in this case S A . The quantity 75' 

denotes the situation where the strobe track only points to the target in question. 
Since the crosses are mutually exclusive, result 2 of Appendix A can be used in 
order to obtain 



P(S^TX)=P(SJ7X;V...VTX„ v TS) 




M 



J^P{TX k ) + P{TS) 



which gives the end result 




Y,P{TX\) + P(TS) 
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^Here, the quantities p[s A \TXl) and p{sjTs),.k = 1 M can be computed exactly 

as in (8). The quantities P(TX' k ), k = i,...,A/ are the prior probabilities of a true 

cross. These are to be selected by the user keeping in mind that the prior 
probabilities of a true cross and P(TS) shall sum up to one along each strobe track 
(including the cross under evaluation). Examples of possibilities are to set the all 
equal or equal to P(order(X)) i.e. to a quantity that depend on the order of the 

cross. In this way higher order crosses could be further favored. As an other 
alternative, a dependence on the geographical location of the cross could be 
introduced. The quantity P{TS) is the prior probability that the strobe track 
corresponds only to a direction. Also this quantity is at the users disposal. 

All procedures necessary to update the quality of a cross X to include the 
information of the strobe track .V is now available. At the beginning the quantity 

p[jX\S h iSvJ is available from previous iterations, performed exactly as the last 

one. The target type probabilities of the cross under evaluation is available from the 
calculation of crosses described above. The target type probabilities of the crosses, 
in which the strobe track .y participates, are also available. The prior probabilities 

p{TX\) and P{TS) are also available. 

The inclusion of the information of the strobe track is then as follows. First 
p[s h \TX,S JiS ...,S Jtt ) is computed according to (8). Then P^^TX'^k = 1 M are 

computed using the technique of (7). Then p[sj-JX y S j2 S Jh ) can be evaluated 

using the prior probabilities. All quantities are now available for the computation of 
piyTX\S h ,...,S >w ), said computation being performed according to (6). 
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So far only the inclusion of the final strobe track in the cross quality has been 
described. Since the inclusion is formulated recursively, a recursive computation 
scheme can now be formulated for the full evaluation of the cross quality. It is 
assumed that all crosses has been computed as described above. It will be obvious 
for anyone skilled in the art that the recursive cross quality computation can be 
formulated as follows; 

1. (Prerequisites): Select prior cross probabilities, P{TX) t for all existing crosses 
and P{TS) for all existing strobe tracks. 

2. (Evaluation of the quality of cross X, P Q [TX\S h S Jk )y 

P Q {TX\0) = P(TX) (0 denotes the empty set) 

for m = l,...,/v* repeat ( N is the number of participating strobe tracks) 

p(5 J 7X,^,,..,i- i ) = ip(^|^)/'(^) 

id 

for k - l,..., M repeat (loop is over crosses that have ^participating) 
end 

tiK\ Tx 'M TX '^ +p ( s ^ TS M TS J .) 

j:P{TX t ) + p(TS ]m ) 
NORMALIZATION = p[sjTX y S lmi S yi )/>(7X|^ ,,...,£,,) 

+ hrA\s ; .., ..... 5^1 - p(tX\S^ ,...,s,)) 
. Ms,\tx,s, sAp(tx\s. Sj) 

r Q\' A ^j„<- ••>**) NORMALIZATION 
end 
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it is .obvious to anyone skilled in the art that the above scheme can be simplified 
further by neglecting different dependencies. 

The cross quality computation can be combined with corresponding kinematic cross 
quality computations, e.g. in the preferred multi-sensor multi-target tracking system 
described above. The combination can then be further used in the procedure for 
automatic track initiation in that system. The combined cross quality computation 
can for example be performed by assuming independence between the kinematics 
and the target type information, which leads to a multiplication of the kinematic cross 
quality and the cross quality described in this invention. This combination is 
expected to give significant improvements in the suppression of ghost targets. This 
effect is illustrated in figs. 10a and 10b. 

Note that the quality measuring quantities computed according to (7) and (8) will be 
high only if at least one target type has the property that both the probability of the 
strobe track and the probability of the strobe track cross is high. Otherwise a low 
quality will result. This is desirable, since if there were no target type with the above 
property, then the computed strobe track cross would be inconsistent with the 
information of this participating strobe track. This effect is illustrated in fig. 9a and b, 
where fig. 9a shows consistency and fig. 9b shows inconsistency. 

In fig. 9a, a strobe track 1 has a high probability for two of the target types, namely 
type 2 and 3. The 2 nd order strobe track cross have only a high probability for the 
target type 2. However, this is enough to produce a high strobe track cross quality, 
since target type 2 also had a high probability in the strobe track 1. In fig. 9b, the 
same strobe track is used. Here, another 2 nd order strobe track cross is considered, 
which have a high probability only for target type 1. In this case, the strobe track 
cross quality will be low, since there is no agreement between the probability 
magnitudes. 
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To further illustrate the operation of the recursive cross quality evaluation discussed 
above an example is given with reference to figs. 10a and 1 0b- 
Example 5 

In fig. 10a a target tracking system without target type information is shown. Without 
target type probabilitiy estimates of the strobe tracker, there is no way to distinguish 
between true and ghost targets in this 2 sensor scenario with azimuthal only strobe 
tracks. All crosses are equally good. In fig. 10b a target tracking system with target 
type information is shown. Here three true targets, denoted by arrow symbols, are 
tracked by 2 sensors S1 and S2 that each sees all targets. As can be seen there are 
8 crosses, denoted by open circles, whose quality needs to be evaluated. These are 
given the numbers 1 to 8 with the numbering running from left to right for each row 
and then from bottom to top. With this numbering the true crosses have numbers 4 
( T x ), 5 ( T 2 ) and 6(7:). Target number 1 thus has true type 7; , target number 2 has 
true type T 2 and target number 3 has true type 7 3 . The strobe tracks of each sensor 
are numbered clockwise starting with sensor 1 (S1). The strobe tracks are thus 
numbered from 1 to 6. There are three possible target types in the scenario. The 
starting point is the following target type probabilities for each strobe track; 
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Using the final result above to compute the crosses (all prior target type probabilities 
are chosen equal, i.e, p(t*)= 1/3) gives 
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The quality of the crosses are then evaluated. It is assumed that the following 
values are used 

P{TX) - 0.25 ; P(TS) = 0.25 for strobe tracks with three crosses. 

P{TX) = 0.33 ; P(TS) = 0.33 for strobe tracks with two crosses. 

This resulted in the following values of p(TX\S h when the recursions where run 

X r 0.1649 
Ay. 0.1 649 
X, : 0.1 649 
X A \ 0.4898 
. X % : 0.4059 
X 6 : 0.4997 
X 7 : 0.2039 
^ 8 :0.2169 

For each cross the strobe track with the least number of crosses was included first. 
Clearly, the correct targets are detected by the method according to the invention. 

Even if the above examples all refer to one particular described MST system, the 
target type estimation method according to the present invention may be used in 
other tracking systems as well. The method is applicable in all cases where tracks of 
different types are formed and/or where these tracks are used to create other 
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related quantities. The discrete information may in this sense be integrated in many 
type of situations. The scope of the present invention is solely determined by the 
attached claims. 
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Appendix: Basic results and the calculation of the discrete probability function 

Result 1: Let the events A and B be mutually exclusive and let z be a 
measurement. Then 



DtlA m P(z\A)P(A) + P(z\B)P(B) 
P WW = P{A) + P(B) 



P(:.(AvB)) P(z.Avz.B) 



P(z,A) + P(z.B)- PjzM^B) P(z,A) + P(z.B) P(z\A)P(A) + P(:\B)P(B) 
P(A)+P(B)-P(A,B) ~ P(A) + P(B) " P(A)+P{B) 



This results can be easily generalized to cover an arbitrary number of mutually 
exclusive events. The result is then 

Result 2: Let the events {a„)' x ( be mutually exclusive and let z be a measurement. 
Then 

t,Pl:\A m )P(A m ) 

P(z\A t v/J,v...v/y = ^— ; 

T P (AJ 

The above results are useful in the calculation of cross qualities. 

The case where the two latest measurement outcomes are retained in the target 
tracking process is now discussed. The discrete probability function then becomes 
(it is only known that the target type is one of T, i = I N T .) 
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LV,.pj>) = r{t, r z p (t,)\{T:v...vT;, r ).Z"{t il _ x )) 

/, (^--"( / ,)-{^" v - v7 ^)- zn <v->) 

/^„.{7;"v...v/;; r }.z"( Vl )) 

= / '( / ,--''(',)l{7;"v...v7.;},-"(V l )^"(V:)) 

4,,{y;"v...vr; r },z"(^. ; )) 
x p(, ir {T v...v7;;},z"(v,))' 

('„..) highlights the fact that only previous measurements associated to the track 

in question are considered. By the same motivation as earlier the quotient of the last 
factor in the last equation should tend to 1 as q increases. Hence, using the fact 
that the target types are a set of mutually exclusive and exhaustive gives 



P (<<<:"«A T > V...7;" r }.r"(/,,-.).ZX-:)) 

So 0,.A ») * 1 1 \ x 

/>(/,477v...7:J.rX.,). Z "(V:)) 

^.,--'^,-)-{ 7 i" v -^}- Z " (/ --->) 

/ J (v{7Tv...r; r },z"(V:)) 



SUBSTITUTE SHEET (RULE 26) 



WO 98/58274,.. . PCT/SE98/00956 

53 

• • / ? (/,.."-(/„),{7;"v...v7:; f },r"(v l ),Z"(V,)) 

^(/, ( .{7;-v...v7;" r }.z"(/ iH )) 
£r(>. l --"«. l )\r," .rx.. ).z"c, , ))/'(/.,..-"(/„, )i '/;,". , ))/'(/,. 7*;iz"(/„-, » 



*Z/</,.^(/jir..--(/ < ,))/'(/ v ..-"(/ 1( ,)ir)/'(/,.riz"(v,)). 

Thus a relatively simple sum results. The measurements are now given by 

--"('., )=<K,'. --"(/,.,) = A/;.,. 
This gives the following final result 
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CLAIMS 

1 . A target type estimation method for estimation of target types in a tracking 
system, using discrete information from measurements, said target type estimation 
method operating in a target type system comprising a number N T of target types 
and a number A> w of measurement outcomes, where each target type has a relation 
to at least one of said measurement outcomes and where each measurement 
outcome has a relation to at least one of said target types, said target type 
estimation method comprises an update step including a time recursive calculation 
of target type probabilities, comprising the steps of: 

calculating probabilities ( P(f t( J,\ZO ( ^)) ) at a certain time (/ (/ ) that a track is 
of said target type ( 7] , / = 1 N T ), using measurements ( Z(/ (/ _, ) ) at a previous time 

(',,): 

calculating the likelihood ( l } (t (r z{t q )|7;,Z(/ (H )) ) for a certain measurement 
(-('„)) given the target type and the measurements up to said previous time; and 

calculating probabilities ( p{i tr T\Z{i H )) ) that a track is of said target type at 
the present time from said probabilities ( /'(/„, T\Z{i ir , )) ) that a track is of said target 

type at a previous time (z^,) and said likelihood ( / J (/, r -(> t; )\T t , Z(/ (H )) ); 
characterised in that said update step further comprises the steps of: 

approximating said likelihood ( P{i^z(i q )\T t ,Z(t irl )) ) from a reduced number of 
parameters; and 

restoring ambiguity in said likelihood approximation in cases where at least 
two target types have relations to a different number of measurement outcomes and 
at least one measurement outcome has a relation to at least two of said target type's. 
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2. . A target type estimation method according to claim 1, characterised 
in that said calculation of said probabilities ( T t \Z(t^ ) that a track is of said 
target type at the present time is based on products of said probabilities 
(P^JWZU ,,.,))) that a track is of said target type at a previous time (/ ir , ) and said 

likelihood (r(i tr =(t it )\r„z(t tH ))); 

3. A target type estimation method according to claim 1 or 2, 
characterised in that said calculation of said probabilities 

( r[t ir 7 ( '|Z(/ (/ )) ) that a track is of said target type at the present time is based on the 
following : 

I^CK-ze.^Mv;! 2 ".-.)) 

1=1 

p{t u j:\7Ai lt ))= r;\ for/ = i * ri 

where P" are the prior target type probabilities. 

4. A target type estimation method according to anyone of the claims 1 to 3, 
characterised in that said approximation comprises the step of 
neglecting all but a fixed number of the previous measurements in the calculation of 
said likelihood. 

5. A target type estimation method according to claim 4, characterised 
in that said approximation comprises the step of neglecting all previous 
measurements except the two last measurements in the calculation of said 
likelihood. 
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6. . A target type estimation method according to any of the preceding claims, 
characterised in that said ambiguity restoring comprises the further 
steps of: 

determining probabilities for a designing target type ( J D )\ and 
using said probabilities obtained for a designing target type ( */* ) as 
probabilities for target types presenting an ambiguity of measurement outcomes. 

7. A target type estimation method according to claim 6, in which said target 
type system is defined by a target type system matrix A rA ' , in which target types are 
running along the columns, the measurement outcomes along the rows and in which 
a relation between a target type and a measurement outcome is defined with 1 and 
a non-existing relation with 0; characterised by in said update step 

approximating said likelihood (/'(/,,,::(/., )|7;.Z(/, , ))) based on said matrix A™ . 

8. A target type estimation method according to claim 7, characterised 
b y in said update step approximating said likelihood ( /'(/„, r(/ v )| 7;, Z(/ v _, )) ) based 
on one single parameter and said matrix A n/ . 

9. A target type estimation method according to claim 8, characterised 

b y in said update step approximating said likelihood (P(f tr z{/ it )\T t ,Z(/ ti _ ] )^) as a 

product of two factors, where the first factor is given by a parameter representing the 
probability of correct detection (P r ) in cases where both the measurements z(t H ) 

and z(t ti , ) are equal and consistent with the target type ( T t ), by a parameter 
representing the probability of erroneous detection ( P E ) in cases where the 
measurement r(/ w .,) is consistent with the target type ( 7j), and the measuremenf " 
r(/ y ) is different from the measurement r(/ iH ), by said parameter ( P k . ) in cases 
where the measurement r(/ , ) is inconsistent with the target type ( T i ) and where 
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the measurement r^) also is inconsistent with the target type ( T ) , and by a 
parameter ( P f ) representing the probability of a consistent detection for the 
designing target type ( 7^) in cases where the measurement is inconsistent 

with the target type ( T. ) and where the measurement z(/ i} ) is consistent with the 
target type ( T t ) t and where the second factor is given by said parameter ( I\ ) when 
the measurement r(/ v .,) is consistent with the target type ( T ), by said parameter 
( I), ) when the measurement r(/^,) is inconsistent with the target type ( 7j ), whereby 

I - /> 

I> = — 

N At - I 

and 

where N mU (tf) denotes the number of measurement outcomes for said designing 
target whereby one of said parameters ( /{. , I), , ) constitutes said single 
parameter. 

1 0. A target type estimation method according to claim 9, characterised 
b y in said update step approximating said likelihood by the product : 

)| 7,\Z( )) = A/J'/;, A/,.,)/^/,, A/^l '/;) , / =l N T 

where the measurement ) is given by the measurement outcome {M ) 
and where the measurement z(t ,) is given by the measurement outcome 
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. and where the factors of said product are given by: 
for/=l M^=l N„ 

where the subscripts and u denotes the row /, and column /. , respectively, 

of a matrix, 
and where 

l\<^\T, , A*,.,) = (J>. - /> K v , , • A™ • A™ ,+!>■ A™ (l - A~ .,) + 

for / = l N T% M H = 1 N At . M q , = ! N K} and where J ( is a 

Kronecker delta function. 

11. A target type estimation method according to any of the preceding claims, 
characterised by using detected ESM emitter modes as measurements. 

12. A target type estimation method according to any of the claims 1 to 10, 
characterised by using detected ESM emitters as measurements. 

13. A target type estimation method according to any of the claims 1 to 10, 
characterised by using detected ESM target types as measurements. 

14. A target type estimation method according to any of the claims 1 to 10, 
characterised by using IRST information as measurements. 

1 5. A target type estimation method according to any of the preceding claims, 
characterised by using outcomes from more than one type of 
measurements, whereby the different target typeprobabilities are based on at least 
two of the available types of measurements. 
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16. A target type estimation method according to claim 1 5, 
characterised by using measurements of at least one of ESM emitter 
modes, ESM emitters, ESM target types, IRST information and direct observations. 

17. A target type estimation method according to any of the preceding claims, 
characterised in that said method further comprises the step of time 
propagation of target type probabilities to the time of the next measurement. 

18. A target type estimation method according to claim 17, 
characterised in that said time propagation of target type probabilities 
is performed according to: 

for / = ! N Ti 

where r is the decay time. 

19. A target type estimation method according to any of the preceding claims, 
characterised by using said target type estimation method in a strobe 
track. 

20. A target type estimation method according to any of the claims 1 to 18, 
characterised by using said target type estimation method in a system 
track. 

21 . A target tracking method, characterised in that it comprises the 
step of calculating the goodness of fit between discrete target type related 
measurements and related estimated target type probabilities, said target 
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.tracking method operating in a target type system comprising a number A'* r of target 
types and a number A' A/ of measurement outcomes, where each target type has a 
relation to at least one of said measurement outcomes and where each 
measurement outcome has a relation to at least one of said target types, said step 
of calculating the goodness of fit comprises the step of calculating a discrete 
probability measurement function (f l} (t r pji)) t where /; refers a particular 

measurement and // to a particular track, where 



said step of calculating said discrete probability measurement function (/ 0 U^p,n)) 

comprises the step of approximating said discrete probability measurement function 
(/„(/„,/>.«)) from a reduced number of measurements. 

22. A target tracking method according to claim 21, characterised in 
that said approximation comprises the step of neglecting all but a fixed number of 
the previous measurements in the calculation of said discrete probability 
measurement function (./*„(/,,,/>,»)) 

23. A target tracking method according to claim 22, characterised in 
that said approximation comprises the step of neglecting all previous 
measurements except the two last measurements in the calculation of said discrete 
probability measurement function (f u (t tr pji))- 

24. A target tracking method according to claim 23, characterised in_ 
that said step of calculating the goodness of fit comprises the step of 
approximating said discrete probability measurement function (f u (t r pji)) as a sum 
of products of three factors, according to 
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A' - 

/ D h-A'.»)=i^.w;ir.w;,M/,.^; 1 ir)/ , (',.riz"(/,-,)). 

where the first factor is given by a parameter representing the probability of correct 
detection ( P c ) in cases where both the measurement outcomes ( M* , M n x ) are 
equal and consistent with the target type ( T ), by a parameter representing the 
probability of erroneous detection ( P } , ) in cases where the measurement outcome 
A-/" , is consistent with the target type ( 7\ ), and the measurement outcome Mf is 
different from the measurement outcome A/" , t by said parameter ( P k . ) in cases 
where the measurement outcome M" t , , is inconsistent with the target type ( T ) and 
where the measurement outcome also is inconsistent with the target type ( T ) , 
and by a parameter ( I) ) representing the probability of a consistent detection for the 
designing target type ( t d ) in cases where the measurement outcome is 
inconsistent with the target type ( '/* ) and where the measurement outcome Mf is 
consistent with the target type ( 7] ) t where the second factor is given by said 
parameter ( P t } when the measurement outcome A-/J , is consistent with the target 
type ( '/' ), by said parameter ( P k . ) when the measurement outcome M n M , x is 
inconsistent with the target type ( 7\ ), and where the third factor is the estimated 
target type probability of target type ( 7] ) of track , whereby 

p 

' : N u -\ 

and 
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A' A ,(</)(* A , -1) 

where N Si (J) denotes the number of measurement outcomes for said designing 
target r u , whereby one of said parameters (P c , I\. , />, ) is selectable by the user. 

25. A target tracking method according to any of the claims 21 to 24, 
characterised in that it further comprises the step of associating a 
discrete measurement (r(/ (/ )) to a track, comprising the step of calculating the 

goodness of fit between said discrete measurement (z(f tl )) and related estimated 
target type probabilities of said track, based on the results from the step of 
calculating a discrete probability measurement function (//,(/<,,/>,«))■ 

26. A target tracking method according to claim 25, characterised in 
that said associating step comprises at least one of the steps of: 

extending each track with a discrete state comprising target type 
probabilities; 

updating and propagating each track to the next measurement time; 
computing a discrete probability function (f^t^pji) ) for any measurement to 

track association; and 

calculating an association quality quantity for each measurement to track 
association based on said discrete probability function (f u (f q ,pjf)). 

27. A target tracking method according to claim 26, characterised in 
that said associating step comprises the step of: ' 

calculating an association quality quantity for each measurement to track 
association based on both said discrete probability function (/ o 0 tr £>,//)) and a 

kinematic only probability density function (A(/ (/ , 
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28. A target tracking method according to claim 26 or 27, characterised 
in that said associating step comprises the step of: 

calculating an association quality quantity for each measurement to track 
association comprising a logarithmic probability quotient (LPQ). 



29. A target tracking method according to claim 28, characterised in 
that said logarithmic probability quotient (LPQ) is calculated according to: 



LPO = 



where 



tp<f(f , P. ft) = lorn 

' ' " ' \ p s (\- p^pcri)) 
where 

P t) is the probability of detection of the sensor, P s is the probability of false 
alarm of the sensor and P(T! ) is the probability that the track corresponds to a true 
target. 



30. A target tracking method according to any of the claims 21 to 29, 
characterised by the step of associating a discrete measurement (z(t q )) 

to a strobe track. 



31 . A target tracking method according to any of the claims 21 to 29, 
characterised by the step of associating a discrete measurement ( z{/ tj ) \ 

to a system track. 



SUBSTITUTE SHEET (RULE 26) 



WO 98/58274" PCT/SE98/00956 



64 



32. - A target tracking method according to any of the claims 21 to 31, 
characterised by the step of calculating a track quality. 

33. A target tracking method according to claim 32, characterised by 
the step of calculating a track quality using said discrete probability function 

34. A target tracking method according to claim 33, characterised by 
the step of calculating a track quality using both said discrete probability function 
(./„(/,,, />,")) and a kinematic only probability density function if K (f tr pji)) 

35. A target tracking method according to claim 34, characterised by 
the step of calculating a track quality according to : 

WV,/ ) = 



i + _ P s l-/ 3 (77\y,) 

^;A(V p. " )f t , C„ ,/>.//)+o- p u )/;. /'( Tl\ / lH ) 

when a measurement is obtained during the scan and according to: 



' !-''„ «7T./,,) 



when no measurement is obtained durina the scan. 



36. A target tracking method according to any of the claims 32 to 35, 
characterised by the step of evaluating said track quality in multiple 
hypothesis tracking. 
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^37. A target tracking method according to any of the preceding claims, 
characterised in that target type probabilities are calculated according 
to a target type estimation method for estimation of target types in a tracking system, 
using discrete information from measurements, said target type estimation method 
operating in a target type system comprising a number N r of target types and a 
number A' w of measurement outcomes, where each target type has a relation to at 
least one of said measurement outcomes and where each measurement outcome 
has a relation to at least one of said target types, said target type estimation method 
comprises an update step including a time recursive calculation of target type 
probabilities, comprising the steps of: 

calculating probabilities ( p{t tr T,\Z(/ tl , )) ) at a certain time ( ) that a track is 

of said target type ( '/; , / = ! N 1 ), using measurements ( Z(/ v , ) ) at a previous time 

(',,): 

calculating the likelihood ( p[t ir r(/ (/ )]'/;, Z(/ (/ , )) ) for a certain measurement 
(-('.,)) given the target type and the measurements up to said previous time; and 

calculating probabilities ( p[t ir T f \Z(/ t) )) ) that a track is of said target type at 

the present time from said probabilities ( p[i ti . T\Z(i u , )) ) that a track is of said target 

type at a previous time (/ (/ t ) and said likelihood ( r(/ v )\T t ,Z(f it _ ] )) ); 

approximating said likelihood ( r(/ (/ )|7,\Z(/ (H )) ) from a reduced number of 
parameters; and 

restoring ambiguity in said likelihood approximation in cases where at least 
two target types have relations to a different number of measurement outcomes and 
at least one measurement outcome has a relation to at least two of said target types. 

38. A cross calculation method, characterised in that it comprises 
the step of calculating the target type probabilities ( P{jl' v \S h A\ v ) , / - I A ; 7 ) of 
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3. cross ( X ) consisting of a subset of alt available strobe tracks (.V, S ), the 

number of strobe track being at least two : according to: 



where the quantities P\S h \ T t x ) , k = 1 /V , are the target type probabilities of the 

strobe track and where / ) (/; v ), k = 1 A' , are the prior target type probabilities of 

the cross. 

39. A cross calculation method according to claim 38, characterised by 
the step of calculating a discrete quality of the cross ( X ), said discrete quality being 
the probability of a true cross conditioned on all participating strobe tracks 

( p['IX\S h ,V ;v ) ) used in the computation of said cross. 



40. A cross calculation method according to claim 39, characterised in 
that the inclusion of an additional strobe track in said discrete quality is 
performed according to: 



''K s ,. s tK ) 

= 4vj7y,.v ft ....,.s; t )^7y|.v J; .... 1 .v te ) 




and where the following approximations are used: 
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41 . A cross calculation method according to claim 40, characterised in 
that the calculation of said discrete quality is performed according to the following 
steps: 

1. selecting prior cross probabilities ( P{7X) ) for all existing crosses and prior strobe 
track probabilities (i } (TS)) for all existing strobe tracks, said prior cross 
probabilities and said prior strobe track probabilities summing up to one along 
each strobe track; 

2. evaluating the quality of cross (A'), {P (J {TX\S h S\ J): 

calculating P Q {'lX\0) = P{TX) , where 0 denotes the empty set; 
repeating the steps of: 

calculating p(s\ m | TX.S, m _ S lt ) = £ P(s im 1 7,' v )p{ T* ) 

calculating p(s ' u \ TX t ) = £ l'(s u \ if' )p(lf ) for k = I M 

calcu lating p(s u | -, TX, S.^ S h ) 

T i^sjixjiiw) * p(s jrs J/'(?:v J 

_ j_t 

X/'(7X) + /'(7-.vJ 

calculating NORMA LIZA TION = p{s im \ TX,S, m ( .V ;i )p(rX\S lm , S J 

*P(sj^TX.S JmJt S^-lipM^ Sj) 
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" , x Ms. \ rx,s l s j )p(tx\s j s.) 

ca,cu,a,in g , c , ( W |.V,. s„). ' 

for m = 1 N , 

where P{lX\ ) is a prior probability of a true cross. 

42 A system for target tracking, characterised by a cross calculation 
step, comprising the step of : 

calculating the target type probabilities (P(7' l x \S h S Jk ) t / = lr....A' 7 ) of a 

cross ( X ) consisting of a subset {S A' A ) of all available strobe tracks, the 

number of strobe tracks being at least two, according to: 



, ('t'lrn+J': 1 '] 



where the quantities P(sjT t x ), k = i N , are the target type probabilities of the 

strobe track and where / J ('/' v ) , k - I /V , are the prior target type probabilities of 

the cross; 

calculating of a discrete quality (P (J [lX\S Jt .Y ; J) of a cross ( A') f said 

discrete quality being the probability of a true cross conditioned on the strobe tracks 
that are participating in said cross ( X ); 

calculation of said discrete quality is performed according to the following 

steps: 

1 . selecting prior cross probabilities ( P{TX) ) for all existing crosses and prior strobe 
track probabilities (P(TS)) for all existing strobe tracks, said prior cross 



SUBSTITUTE SHEET (RULE 26) 



WO 98/58274, , 



PCT/SE98/00956 



69 

- probabilities and said prior strobe track probabilities summing up to one along 
each strobe track; 

2. evaluating the quality of cross ( X ), (P q [tx\S u S Js )): 

calculating P q {'IX\q) = P{TX) , where 0 denotes the empty set; 
repeating the steps of: 

calculating />(.Y,_ ! TX. S Jm _ S J = £ 4 V i- 1 7 .' v H^ 1 ') 

calculating p(s jm | 'A', ) = £ 7 i V y. I 'C v; M'A""' ) for * = I A/ 

calculating p[s. m j -. 7*, ^ ( A" J 

tp(rx l )^p('rs jm ) 

calculating NORMA LIZA TION = l>(s im \ TX. S Jm t V t )p(rX\S Ui 5 J 

+ p(sj-.TX.s\ mi s h )(\-p(rx\s^ i'J) 

calculating pJ/XlS S )= 1 ; " . im> hj V ^ 

for /77 = I yV , 

where P{TX\) is a prior probability of a true cross; 

and a step of estimation of target types in a tracking system, using discrete 
information from measurements, said target type estimation step operating in a 
target type system comprising a number A', of target types and a number A' A/ of 
measurement outcomes, where each target type has a relation to at least one of 
said measurement outcomes and where each measurement outcome has a relationr 
to at least one of said target types, said target type estimation step comprises an 
update step including a time recursive calculation of target type probabilities, 
comprising the steps of: 
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- calculating probabilities ( P(t q J,\Z(t iH )) ) at a certain time ( / (/ ) that a track is 
of said target type ( T , / = I N T ), using measurements ( at a previous time 

(',.,); 

calculating the likelihood ( /'(/„.-(/„ )\T t , Z(t H . x ))) for a certain measurement 
(r(/ t; )) given the target type and the measurements up to said previous time; 

calculating probabilities {P(f q ,T l \Z(f ti )) ) that a track is of said target type at 

the present time from said probabilities ( p[t^ 7;|Z(/ (/ _, )) ) that a track is of said target 

type at a previous time and said likelihood ( )| 7), Z(/ v ., )) ); 

approximating said likelihood ( ^(/^(/JIT^Z^.,))) from a reduced number of 
parameters; and 

restoring ambiguity in said likelihood approximation in cases where at least 
two target types have relations to a different number of measurement outcomes and 
at least one measurement outcome has a relation to at least two of said target types. 

43 A system for target tracking according to claim 42, characterised 
b y a step of calculating the goodness of fit between discrete target type related 
measurements (z(r q )) and related estimated target type probabilities, said step of 

calculating the goodness of fit comprises the step of calculating a discrete 
probability measurement function (f t) V r pjt)), where /; refers a particular 

measurement and // to a particular track, where 

said step of calculating said discrete probability measurement function (/ D (/„,/>.»)) 
comprises the step of approximating said discrete probability measurement function 
(//;('. r / ; '")) fro™ a reduced number of measurements. 
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